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POPULAR SUMMARY 


Variations in Stratospheric inorganic Chlorine Between 1991 and 2006 

D.J. Lary, D.W. Waugh, A.R. Douglass, R.S. Stolarski, P.A. Newman, H, Mussa 

So how quickly will the ozone hole recover? This depends on how quickly the chlorine 
content (Cl y ) of the atmosphere will decline. The ozone hole forms over the Antarctic 
each southern spring (September and October). The extremely small ozone amounts in 
the ozone hole are there because of chemical reactions of ozone with chlorine. This 
chlorine originates largely from industrially produced chlorofluorocarbon (CFC) 
compounds. An international agreement, the Montreal Protocol, is drastically reducing 
the amount of chlorine-containing compounds that we are releasing into the atmosphere. 

To be able to attribute changes in stratospheric ozone to changes in chlorine we need to 
know the distribution of atmospheric chlorine. However, due to a lack of continuous 
observations of all the key chlorine gases, producing a continuous time series of 
stratospheric chlorine has not been achieved to date. We have for the first time devised a 
technique to make a 17-year time series for stratospheric chlorine that uses the long time 
series of HC1 observations made from several space borne instruments and a neural 
network. The neural networks allow us to both inter-calibrate the various HC1 
instruments and to infer the total amount of atmospheric chlorine from HC1. These new 
estimates of Cl y provide a much needed critical test for current global models that 
currently predict significant differences in both Cl y and ozone recovery. These models 
exhibit differences in their projection of the recovery time and our chlorine content time 
series will help separate the good from the bad in these projections. 
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3 A consistent time series of stratospheric inorganic chlorine Cl y from 1991 
i to present is formed using space-borne observations together with neural net- 
s works. A neural network is first used to account for inter-instrument biasses 
e in HC 1 observations. A second neural network is used to learn the abundance 
7 of Cl y as a function of HC 1 and CH4, and to form a time series using avail- 
s able HC 1 and CH 4 measurements. The estimates of Cl y are broadly consis- 

9 tent with calculations based on tracer fractional releases and previous esti- 

10 mates of stratospheric age of air. These new estimates of Cl y provide a crit- 
n ical test for current global models that predict significant differences in Cl y 
12 and ozone recovery. 
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1. Introduction 

is Knowledge of the distribution of inorganic chlorine Cl y in the stratosphere is needed to 
H attribute changes in stratospheric ozone to changes in halogens, and to assess the realism 
is of chemistry-climate models [ Eyring et al, 2006; Eyring, 2007], However, there are limited 
is direct observations of Cl y . Simultaneous measurements of the major inorganic chlorine 
1 ? species are rare [Zander et al, 1992; Gunson et al, 1994; Bonne et al, 2000; Nassar et al, 
is 2006]. In the upper stratosphere, Cl y can be inferred from HC1 alone (e.g., Anderson et al 

i, [ 2000 ]). 

20 Here we combine observations from several space-borne instruments using neural net- 

21 works [ Lary and Mussa, 2004] to produce a time series for Cl y . A neural network is used 

22 to characterize differences among various HC1 measurements, and to perform an inter- 

23 instrument bias correction. Measurements from several different instruments are used in 

24 this analysis. These instruments, together with temporal coverage and measurement un- 

25 certainties, are listed in Table 1. All instruments provide measurements through the depth 

26 of the stratosphere. A second neural network is used to infer Cl y from these corrected 

27 HC1 measurements and measurements of CH 4 . 

28 Sections 2 and 3 describe the HC1 and Cl y inter comparisons. Section 4 present a sum- 

29 mary. 

2. HC1 Intercomparison 

30 We first compare measurements of HC1 from different instruments listed in Table 1. 

31 Comparisons are made in equivalent PV latitude - potential temperature coordinates 

32 [Schoeberl et al, 1989; Proffitt et al, 1989; Lait et al, 1990; Douglass et al, 1990; Lary' 
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33 et al, 1995; Schoeberl et al, 2000] to extend the effective latitudinal coverage of the 
3 i measurements and identify contemporaneous measurements in similar air masses. 

35 The Halogen Occultation Experiment (HALOE) provides the longest record of space 

36 based HC1 observations. Figure 1 compares HALOE HC1 with HC1 observations from 
a? (a) the Atmospheric Trace Molecule Spectroscopy Experiment (ATMOS), (b) the Atmo- 

38 spheric Chemistry Experiment (ACE) and (c) the Microwave Limb Sounder (MLS). In 

39 these plots each point is the median HC1 observation made by the instrument during each 
lo month for 30 equivalent latitude bins from pole to pole and 25 potential temperature bins 

4 1 from the 300-2500 K potential temperature surfaces. 

42 A consistent picture is seen in these plots: HALOE HC1 measurements are lower than 

43 those from the other instruments. The slopes of the linear fits (relative scaling) are 
« 1.05 for the HALOE-ATMOS comparison, 1.09 for the HALOE-MLS, and 1.18 for the 
is HALOE-ACE. The offsets are apparent at the 525 K isentropic surface and above. Pre- 
46 vious comparisons among HC1 datasets reveal a similar bias for HALOE [Russell et al, 
i? 1996; McHugh et al., 2005; Froidevaux et al. , 2006]. ACE and MLS HC1 measurements are 

48 in much better agreement [Figure 1(d)]. Note, all measurements agree within the stated 

49 observational uncertainties summarized in Table 1 . 

so To combine the above HC1 measurements to form a continuous time series of HC1 (and 
si then Clj,) from 1991 to 2006 it is necessary to account for the baises between data sets. A 

52 neural network is used to learn the mapping from one set of measurements onto another as 

53 a function of equivalent latitude and potential temperature [ Lary and Mussa , 2004], We 

54 consider two cases. In one case ACE HC1 is taken as the reference and the HALOE and 


DRAFT 


April 11, 2007, 3:27pm 


DRAFT 



LARY ET AL.: STRATOSPHERIC CHLORINE FROM 1991-2006 X - 5 

55 Aura HC1 observations are adjusted to agree with ACE HC1. In the other case HALOE 
ss HC1 is taken as the reference and the Aura and ACE HC1 observations are adjusted to agree 
s? with HALOE HC1. In both cases we use equivalent latitude and potential temperature 
ss to produce average profiles. The purpose of the mapping is simply to learn the bias as a 
59 function of location, not to imply which instrument is correct. 

m The precision of the correction using the neural network mapping is of the order of ss 

6 1 0.3 ppbv, as seen in Figure 1(e) which shows the results when HALOE HC1 measurements 

62 have been mapped into ACE measurements. The mapping has removed the bias between 

63 the measurements and has also straightened out the ‘wiggles’ in 1 (c), i.e. , the neural 

64 network has learned the equivalent PV latitude and potential temperature dependence 
es of the bias between HALOE and MLS. The inter-instrument offsets are not constant in 
66 space or time, and are not a simple function of Cly. 

3. Inorganic Chlorine Cly 

er To a first approximation Cly HC1 + CIONO 2 + CIO [Brasseur and Solomon, 1987], 
es and Cl y can be estimated from HC1 and CIONO 2 . However, observations of CIONO 2 are 

69 much more limited than from HC1. As shown in Table 1, CIONO 2 measurements have 

70 been made by the CLAES (1991-1993), ATMOS (1992-1994), CRISTA (1994, 1998), and 

71 ACE (2004-present). 

72 Because of the limited temporal coverage of CIONO 2 measurements it is not possible 

73 to form a continuous time series of Cl y by combining HC1, CIONO 2 , and CIO. However, 

74 it is possible to form a time series of Cl y using a neural network. There are sufficient 
7s observations of CIONO 2 from ATMOS, CLAES, CRISTA, and ACE to train a neural 
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76 network to learn the Cl y abundance as a function of HC 1 and CH 4) for each of which there 

77 is a long, near-continuous, time series of measurements. The resulting reconstruction 

78 reproduces an independent validation dataset faithfully with a correlation coefficient of 

79 0 . 99 , and provides a scatter diagram with a slope very close to one for the observed Cly 
so plotted against the neural network inferred Cl y , see Figure 1(f). 

si The inputs to the neural network that estimates Cl y are HC 1 , CH4, equivalent latitude 
82 and potential temperature. HC 1 is used because it is continuously observed from the 
as launch of UARS to the present and is typically the major Cl y reservoir. CH4 is used 
84 because it is continuously observed from the launch of UARS to the present and, as a 
ss long-lived tracer, it is well correlated with Cl y . Potential temperature and equivalent 
as latitude are used because the correlation between long-lived tracers such as CH 4 and Cl y 
s? is a strong function of altitude and a weak function of latitude [ Lary and Mussa, 2004 ]. 
ss Other training strategies using more species were examined. For example, we tested 
as the effectiveness of a neural network with inputs of HC 1 , O3, CH 4 , H2O , equivalent 

90 latitude and potential temperature to estimate Cl y . This was tried as O3, CH 4 and H2O 

91 are key observed species involved in the partitioning of reactive chlorine. When chlorine 

92 atoms are released from the chlorine containing source gases by photolysis, they react 

93 with CH 4 to form HC 1 . Alternatively, Cl atoms may react with ozone to form CIO, and 

94 then CIO will combine with NO2 to form C 10 N 0 2 . HC 1 is destroyed either by reaction 

95 with OH, photolysis or heterogeneous reactions. The amount of OH present depends on 

96 the photolysis of ozone to form 0 ( X D) and the subsequent reaction of 0 ( 1 D) with H2O. 
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97 This approach also gave good results, but with slightly lower skill than just using HC1, 

98 CH 4 , equivalent latitude and potential temperature to estimate Cly. 

99 Figure 2 shows how Cl y profiles estimated by the neural network agree with observed 

100 Cly for October 2006. In each case the shaded range represents the uncertainty associated 

101 with the Cl y estimate. We note that the HC1 bias between HALOE and ACE is the major 

102 uncertainty. 

103 The distribution of Cl y is expected to change between 1991 and 2006 as the abundances 

104 of its source gases have changed. Figure 3 shows the time-series of Cly for the 525 K 

105 isentropic surface (f» 20 km) and the 800 K isentropic surface (rs 30 km), for three 

106 different equivalent latitudes. The upper limit of each shaded range corresponds to the 

107 estimate of Cly for the neural network calibrated to agree with ACE v2.2 HC1, and the 

108 lower limit to the estimate of Cly for the neural network calibrated to agree with HALOE 

109 vl9 HC1. 

no The variation in Cly estimates between the two cases depends on latitude, altitude 
m and season and is typically <0.4 ppbv at 800 K. This uncertainty is primarily due to 
in the discrepancy between the different observations of HC1 which translates into the Cly 
113 uncertainty shown by the shading in Figure 3. There is also a slight low bias in the lower 
in stratosphere due to not including HOC1 in the estimates of Cly. HOC1 was not included 
ns because HOC1 has been observed by ACE only since the start of 2004. Ignoring HOC1 
ns is only of significance in regions of strong chlorine activation at low temperatures in the 
a? lower stratosphere where HOC1 can comprise up-to about 10% of Cly. 
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ns There is a general tendency of Cl y to increase in the 1990s, peak around 2000, and 
ns then slowly decrease. This is consistent with our expectations based on the tropospheric 

120 abundence of chlorine containing source gases. The Cl y time-series shown in Figure 3 

121 constitutes a useful test for model simulations. The variation in simulated Cl y from the 

122 chemistry-climate models used in the recent WMO [2006] report is much greater than 

123 the above uncertainty in Cl y . For example, the simulated peak Cl y in October at 80S 

124 varies from less than 1 ppbv to over 3.5 ppbv, while the peak annual-mean Cl y for north 

125 mid-latitudes varies from 0.8 to 2.8 ppb [ Eyring et al, 2006; Eyring, 2007]. 

we The estimates of Cl y produced are broadly consistent with calculations based on tracer 
127 fractional releases [Newman et al., 2006] and previous estimates of stratospheric age of 
us air. Observations show that at 20 km the mean age increases from around 2 years in 
i29 the tropics to around 4 years at high latitudes (60°N), with a similar gradient at 30 km 
uo but older ages by around 2 years [ Waugh and Hall, 2002], The curves in Figure 3 show 

131 calculations of Cl y for a range values of the mean age of air, and the ages that are required 

132 to match the observed Cl y are consistent with the observations of the mean age. 

4. Summary 

133 A consistent time series of stratospheric Cl y from 1991 to present has been formed 

134 using available space-borne observations. Here we used neural networks to inter-calibrate 

135 HC1 measurements from different instruments, and to estimate Cl y from observations 

136 of HC1 and CH4. These estimates of Cl y peaked in the late 1990s and have begun to 

137 decline as expected from tropospheric measurements of source gases and troposphere to 

138 stratosphere transport times. Furthermore, the estimates of Cl y produced are consistent 
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139 with calculations based on tracer fractional releases and age of air [Newman et al, 2006 ]. 
no The Cl y time-series formed here is an important benchmark for models being used to 
Mi simulate the recovery of the ozone hole. Although there is uncertainty in the estimates 
no of Cly, primarily due to biases in HC 1 measurements, this uncertainly is small compared 

143 with the range of model predictions shown in the recent WMO [ 2006 ] report. The two 

144 Cly time-series are available in the electronic supplement. 

ns Acknowledgments. It is a pleasure to acknowledge NASA for research funding, Lu- 
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Instrument 

Temporal Coverage 

Species 

References 

Median Observation Uncertainty 

ACE 

2004-2006 

HC1, ClONOs and CIO 

Bemath et al. [2005] 

8% (HC1), 30% (C10N0 2 ), >100% (CIO) 

ATMOS 

1991, 1993, 1994 

HC1, C10N0 2 

Zander et al. [1992] 

8% (HC1), 60% (CIONO 2 ) 

Aura MLS 

2004-2006 

HC1 and CIO 

Froidevaux et al. [2006] 

12% (HC1), 76% (CIO) 

CLAES 

1991-1993 

CIONO 2 

Roche et al. [1993] 

>100% 

CRISTA 

1994, 1997 

CIONO 2 

Offermann et al. [1999] 

61% 

HALOE 

1991-2005 

HC1 

Russell et al. [1993] 

4% 


Table 1. The instruments and constituents used in constructing the Clj, record from 
1991-2006. The uncertainties given are the median uncertainties of the level 2 product for 
all the observations made. 
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Figure 1 . Panels (a) to (d) show scatter plots of all contemporaneous observations 
of HC1 made by HALOE, ATMOS, ACE and MLS Aura. In panels (a) to (c) HALOE 
is shown on the x-axis. Panel (e) correspond to panel (c) except that it uses the neural 
network ‘adjusted’ HALOE HC1 values. Panel (f) shows the validation scatter diagram 
of the neural network estimate of Clj, versus the actual Cl y for a totally independent data 
sample not used in training the neural network. 
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1 0/2006 Mean Cly from 30° to 60° 



Figure 2. Cly average profiles between 30° and 60°N for October 2006. The blue 
line shows the Cly estimated by a neural network using HC1 observations calibrated to 
agree with HALOE vl9 HC1. The green line shows the Cl y estimated by a neural network 
using HC1 observations calibrated to agree with HALOE vl9 HC1. The red line shows 
observed Cly=HCl+C10N02+C10 based on ACE v2.2 data. In each case the shaded 
range represents the uncertainty associated with the Cl y estimate. 
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(a) 


(b) 



(c) 

Figure 3. Panels (a) to (c) show October Cl y time-series for the 525 K isentropic 
surface (ps 20 km) and the 800 K isentropic surface (« 30 km). In each case a shaded 
range representing the uncertainty in our estimate of C\ y is shown. This uncertainty is 
due to the biases between the various HC1 observations. The upper limit of the shaded 
range corresponds to the estimate of Cl y based on all the HC1 observations calibrated 
by a neural network to agree with ACE v2.2 HC1. The lower limit of the shaded range 
corresponds to the estimate of Cl v based on all the HC1 observations calibrated to agree 
with HALOE vl9 HC1. Overlaid are lines showing the Clj, based on age of air calculations 
[Newman et ai, 2006]. To minimize variations due to differing data coverage Months 
with less than 100 observations of HC1 in the equivalent latitude bin were left out of the 


time-series. 
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